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^jQ' Abstract. It is well known that, in Landau gauge, the renormalization function of the ghost-ghost- 

^ , gluon vertex Zi (p^) is finite and constant, at least to all orders of perturbation theory. On the other 

• hand, a direct non-perturbative verification of this result using numerical simulations of lattice QCD 

I is still missing. Here we present a preliminary numerical study of the ghost-ghost-gluon vertex 

. and of its corresponding renormalization function using Monte Carlo simulations in SU (2) lattice 

Landau gauge. Data were obtained in 4 dimensions for lattice couplings j3 = 2.2, 2.3, 2.4 and lattice 
sides A? = 4, 8, 16. 

g : INTRODUCTION 

oo . 

O . Although the Faddeev-Popov ghosts are absent from the physical spectrum of QCD, 

^ ! one can use the ghost-ghost-gluon vertex and the ghost propagator to calculate physical 

observables such as the strong running coupling ^ 



„V) = ao?J^S22M. (1) 

Here (Xq = g^/^Tt is the bare coupling constant and Z-i{p^), Z2,{p^) and Z\{p^) are, 
^ ■ respectively, the gluon, ghost and ghost-ghost-gluon vertex renormalization functions. 

■ In Landau gauge the renormalization function of the ghost-ghost-gluon vertex Z\ [p^] 

is finite and constant, i.e. independent of the renormalization scale /i, at least to all orders 
of perturbation theory [2]. This implies that in this gauge one can consider a definition 
of the running coupling constant (in the momentum subtraction scheme) given by ^ 

= ao/£>(/)G2(/), (2) 

where D{p^) and G{p^) are, respectively, the gluon and ghost propagators. An indirect 
evaluation of the ghost-ghost-gluon renormalization function on the lattice has been 
recently presented in [T], confirming that Zi(p^) is finite in the continuum limit. On the 
other hand, a direct non-perturbative verification of this result using lattice simulations 
is still missing. Let us stress that a direct evaluation of Zi {p^) would allow [0] a study of 
the running coupling constant using eq. (HJ) instead of eq. Q. Such a study may improve 
the precision of the determination of as{p^), since in this case one will not need to 
use a matching rescaling technique iQ] when considering data evaluated at different 



lattice couplings /3. Moreover, a direct study of Z\{p^) is interesting to determine its 
j8 -dependence and possible effects due to the breaking of rotational symmetry. Here we 
present a preliminary numerical study of the Landau-gauge ghost-ghost-gluon vertex 
and of its corresponding renormalization function Z\ (p^) in the pure SU (2) case. 



THE VERTEX ON THE LATTICE 

The 3-point function (in momentum space) for ghost, anti-ghost and (gluon) fields is 

V^'^{k;q,s) = (2;r)4 d\k + q-s) {A"^{k)G'^^\q)) , (3) 

where G*'^(y,z) is the inverse of the Faddeev-Popov matrix and 5^{k+q — s) corresponds 
to momentum conservation. Then, we can define the ghost-ghost-gluon vertex function 
by "amputating" the corresponding 3-point function, i.e. 

, V {A"(k)G^''(q)) 

^"^^1 D{k^)G{q^)G{s^) ' ^ ^ 

where V is the lattice volume and s = k + q. From the weak-coupling expansion in 
SU (Nc) lattice gauge theory one obtains that the ghost-ghost-gluon vertex, at tree level, 
is given by L5J 

{k,q) = ig^fb'^q^co^i^-^^^ , (5) 

where a is the lattice spacing. Let us recall that a generic lattice momentum k has 
components k^ = 2sin(^^(3/2)/a, where k^ = Inn^/N^ and n^, A^^ are, respectively, 
the Matsubara modes and the number of lattice points in the direction. Clearly, 
by taking the formal continuum limit a — » in Q one recovers the continuum tree- 
level result Y^^'^{k,q) = igof^^'^q^ ■ Finally, we obtain that the vertex renormalization 
function in momentum- subtraction scheme is given by 



T^ahc 



z 



where Nc is the number of colors and we consider the asymmetric point with zero 
momentum for the gluon = 0), implying s = q. 



NUMERICAL SIMULATIONS AND RESULTS 

For each lattice volume V = A'^^ = 4^, 8"^, 16^ we have performed simulations for three 
lattice couplings (j8 = 2.2,2.3,2.4). We consider lattice momenta with components 
Pi = Pi = P2, = Pa = P and p = Inn/N, where n gets values l,...,N/2. With this 
choice of momenta we obtain 

Zf^/) = . \ ,^, Up') (7) 
^ pcos(p/2) 
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FIGURE 1. Results for ) and Zj (/5 ) with lattice volume V = 16 . For the lattice spacing in 
physical units we refer to [1]. Error bars were obtained using the bootstrap method with 250 samples. 



where we used the fact that only the imaginary part of the vertex function Y'^'^iO.p) 
contributes to Z^^{p^). Results for the quantity E(p^) defined above and the renormal- 
ization function Zi{p^) are presented in Fig.[T] Clearly, T.(p^) has the same momentum 
dependence as the tree-level vertex [i.e. ~ pcos(p/2)], implying that Zi(p^) is approx- 
imately constant with respect to the momentum p. Also the dependence on /3 seems to 
be very weak. More details will be presented elsewhere 
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